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Abstract 

We describe two complementary ways to show the presence of higher order effects in 
the 1/rriQ expansion for inclusive B decays that have been dubbed 'Intrinsic Charm'. 
Apart from the lessons they can teach us about QCD's nonperturbative dynamics their 
consideration is relevant for precise extractions of \ V c b\- for they complement the estimate 
of the potential impact of 1 /mg contributions. We draw semiquantitative conclusions for 
the expected scale of Weak Annihilation in semileptonic B decays, both for its valence 
and non- valence components. 
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1 Introduction 



The Operator Product Expansion (OPE) has turned out to be a highly powerful tool for 
describing inclusive decays of beauty hadrons. Among the major results and successes are the 
precision determinations of CKM elements using inclusive semileptonic decays: e.g., for \V c b\ a 
relative uncertainty of less than 2% has been achieved based on OPE methods. 

To validate and even improve on such accuracy one has to analyze all potential sources of 
uncertainties. The OPE is widely viewed as yielding an expansion in inverse powers of the 
b quark mass m b for sufficiently inclusive quantities. However, in b — > c transitions a second 
scale is present that has a reasonable claim to be considered heavy, namely the charm quark 
mass with m c 3> Aqcd- It is then natural to integrate out the charm quark as well. This is 
usually done for both quarks at the same scale, which means that the ratio m c /mb is treated 
number of order unity. 

On the other hand, numerically we have m 2 c ~ m 6 AQ CD and hence m h 3> m c 3> Aq CD . This 
suggests an alternative where one integrates out the beauty quark in a first step, while leaving 
the charm quark still a dynamical quark; the latter is subsequently removed in a second step. 
In this procedure one has an intermediate theory still containing charm quarks as dynamical 
entities. Such a description contains dimension-six (or higher) operators of the schematic form 



where b v is the (now) nonrelativistic b quark field and T a Dirac matrix depending on the 
process under consideration; it may contain derivatives yielding higher-dimension operators. 
It has been discussed in some detail in [Tj |2] that these operators match onto a contribution 
proportional to the Darwin term p^ D upon integrating out the charm quark. Yet new subtleties 
emerge at this point. This contribution contains an infrared (IR) sensitive piece proportional 
to lnm(,/m c , which diverges for m c — > 0. It is generated in the two-step procedure by the RG 
running between m& and m c , while it appears simply in the coefficient function when integrating 
out bottom and charm together. In fact, it has been shown in pQ that at higher orders in the 
OPE even inverse powers of m c do appear which have a stronger infrared sensitivity. 

Based on a superficial similarity with various non-valence nonperturbative charm effects 
previously discussed in the literature [3], these infrared sensitive terms have been dubbed "In- 
trinsic Charm" (IC) in inclusive B decays. They were largely treated as a factor contributing to 
the theoretical uncertainty from higher-order effects pQ. In this paper we show how the seem- 
ingly different reasonings given in [TJ [2] are equivalent qualitatively as well as quantitatively 
and perform a systematic and more complete treatment of these terms. It turns out that the 
conventional OPE yields a combined expansion in powers of l/m l c with / > (in fact, 

starting with 1 = 2 at tree level, without extra gluon loop) appearing first for k = 3. This gives 
rise to a quantitative subtlety as well: Since as mentioned above m 2 c ~ m&AQ CD , one has to 
account for the terms l/(m\m 2 c ) in order to complement the full set of l/m^ corrections. 

The remainder of the paper is organized as follows. In the next section we perform the 
standard procedure and integrate out bottom and charm together in a similar way at a scale 
below m c and trace the origin of the intrinsic charm contributions in this approach. In Sect.[3]we 
perform the two-step procedure; i.e., first we integrate out just the bottom quark initially leaving 
charm as a dynamical quark in the theory before integrating it out as well at a lower scale. As 
expected these two ways yield the same results equivalent to what has been discussed in [I] . We 
provide a more consistent and concise derivation of the relation between (generalized) Weak 
Annihilation and the four-quark operators than it was done in Ref. jlj. In Sect, plwe discuss 
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explicitly the dimension-eight operators that yield l/m\xl/m 2 terms at tree level qualitatively 
as well as quantitatively, confirming the results of pQ. We summarize the lessons learnt on 
QCD's non-perturbative dynamics, the implementation of the OPE and on the uncertainties in 
the extraction of V cb in Sect. [6] An estimate of the potential scale of Weak Annihilation (WA) 
for b— >u is given there as well. 

2 The standard OPE for B^X c i9 t revisited 

One starts by considering the doubly differential rate for B-^X C lT>i 
d{ fp )dp2 = V(v- P ) 2 -P 2 9({vpf-p 2 ) 9 (m b -vp) 9 {m 2 +p 2 -2m b vp) L^W^(p) 

where p = m b v — q with q denoting momentum transfer to the lepton pair and v the B meson 
velocity. L^ u is the leptonic tensor with 

Lfjvijp) = -ml[g^ u - Vfj,v u ] + m b [2g^(v ■ p) - v^p v - p^v u ] - [g^p 2 - p^p v \ (2.2) 

and W^ v the hadronic counterpart given by 

2M B W (a ,(jp) = J d A x exp(ip • x){B{v)\b v {x)T u c{x) c(0)T^b v {0)\B{v)) (2.3) 

with b v the 'rephased' (nonrelativistic) b quark field, b v (x) = exp(im b (v-x)) b(x) and = 
7m(! -7s)- 

2.1 General argument 

Here we scrutinize the standard way of setting up the OPE for this inclusive process treating 
both the bottom and the charm quarks as heavy and neglecting a hierarchy of scales between 
the bottom and the charm masses. We remove both from the dynamical degrees of freedom in 
one step at a scale fi below the charm mass. No operators with charm fields can then emerge 
in the OPE, and only operators with nonrelativistic b quarks and their (covariant) derivatives 
appear]]] 

The OPE is constructed by contracting the two charm quarks into (the imaginary part of) 
a propagator, while the beauty quark field is replaced by its expansion in inverse powers of m b . 



Formally expanding the bi-local operator in (2.3) in local operators yields, upon integration, 



the well known Heavy Quark Expansion (HQE). The expansion can be done by external field 
methods particularly economical at tree level, see e.g. [5]. This yields (b v = b v (0)) 

2M B W, u (p) = (B(v)\b v r u ^ + m c )T^b v \B(v))5 + (p 2 ~m 2 c ) 

+ {B{v)\b v T v (j) + m c ){ip){i> + m c )TJ> v \B(v)) 5' + (p 2 - m 2 ) 

+ ^{B{v)\b v T v (j> + m c ){ip)(j> + m c ){ilp)(j> + m c )T,b v \B(v)) 5'[{p 2 - m 2 c ) 
+ ... 

00 1 r i n 



J2^-{B{v)\b v r u {p / + m c 



n=0 



(ifl)(j/> + m c 



r»b v \B(v))5™(p 2 -mi). (2.4) 



1 We do not discuss here operators involving light quarks appearing additionally at 0(a s ) level, their contri- 
bution is small. 
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The main challenge in evaluating higher order contributions lies in identifying the independent 
hadronic parameters controlling the 'string' of matrix elements 



{B{v)\b v {iD^)(iD^) ■ ■ ■ {iD, n )b v \B{v)) (2.5) 

entering at order n + 1. It should be noted that these quantities do not depend on p. 

Contracting the hadronic and leptonic tensors for the pseudoscalar B mesons yields a func- 
tion of (y ■ p) and p 2 of the general form 

oo 

L^W* = V(vp) 2 -P 2 Pn((v-p),p 2 )6 { :\p 2 -m 2 c ) (2.6) 

n=0 

with P n (v ■ p) denoting a polynomial in v ■ p. The analysis so far referred to a fully differential 
distribution in general kinematics. To proceed to the inverse mass expansions we need to 
consider partially integrated probabilities. 

The integration over the variable v ■ p has the limits sjp?- < v ■ p < (ml + p 2 ) / (2m{,), which 
yields terms logarithmic in p 2 from the lower end of integration. Focusing on these logarithms 
we get for 1 = 0,1, ... 



d(v ■ p) y/(v ■ p) 2 -p 2 {v- p) 21 = d {p 2 ) l+1 In + • • • , Q = 4V? V(f+2) (2 ' 7) 

d(y ■ p) ^(v-p) 2 -p 2 (v ■ p) 2l+1 = + • • • 



where the ellipses denote polynomial terms in p 2 , and the coefficients Ci are simple fractions: 
C = l/A, Ci = l/16, C 2 = l/32, C 3 = 5/256 ... 

These logarithms are the source of the IR sensitivity of the coefficient functions to the charm 
mass in this approach. This becomes manifest when they are combined with the derivatives of 
the 5-function in Eq.(2.6): 

(p 2 ) k \np 2 5 {k \p 2 -m 2 c ) = (-l) k k\ \nm 2 c 5(p 2 - m\) + • • • (2.8) 

/I \ n ~ k 

(p 2 ) k \np 2 6 {n) (p 2 -m 2 c ) = k\(n-k-l)\ \^—J 5(p 2 -m 2 c ) forn>k (2.9) 

where k is some integer power and the ellipses point to less singular terms as m c — ► 0. Note 
that at the lower limit of integration in Eq. (2.7) we have p\=p 2 , i.e. p^O. This is the infrared 
regime for the charm quark in the final state. 

The leading IR sensitive terms in the integrated rates arise from only the leading term in 
the leptonic tensor 

L^ = -ml[g^-v tl v v \ (2.10) 

already in the expression for the differential rate. For the terms containing the vector p would 
lead to powers of v- p and p 2 . To obtain non-leading in l/m& terms, we consider the subleading 
terms in the leptonic tensor 

L ^ = m b [2 gfiv (v ■ p) - VflPu - Pli v v ] (2.11) 
L^ = -[g^- VlJ p v ]. (2.12) 
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On the other hand, the n th term in the sum (2.4) for the hadronic tensor contains 



P n oc Y v {j) + m c )7 Ml (^ + m c )7 M2 • • • (j/> + m c )7 /Un (^ + m c )Y^ (2.13) 



which yields upon contraction with the leptonic tensor (2.10) and with the nonperturbative 



matrix elements ( |2.5[ ) a contribution of the form (see (2.6)) 



p n = ( v ■ PY(P 2 y( m ty with i + 2j + 2Z = n + 1 . (2.14) 

Note that due to the purely left-handed structure of the current T, we can have only even powers 



of m c here. Furthermore, if n is even, % necessarily is odd (and vice versa). Hence from (2.7) 
we conclude that IR sensitive terms can appear only if n is odd, which means that the number 
of covariant derivatives in the hadronic matrix element has to be odd as well. Therefore all 
operators which contribute to intrinsic charm have to match onto partonic matrix elements with 
at least one gluon. In turn, it implies that there is no intrinsic charm contribution to operators 
of the form {{k 2 ) n ) where k are the spatial components of the residual b quark momentum. 

Now we can trace how singular terms actually emerge. Putting everything together, we get 
for i = 2m 

J d{v ■ p)y/(v ■ p) 2 ~ P 2 Pn{(v ■ p),p 2 )6 (n) {p 2 - m 2 ) (2.15) 
= X)o iil J d(v ■ p)y/(v ■ p) 2 - p 2 [(v ■ P y(p 2 y(m 2 c ) 1 } 5^(p 2 - m 2 c ) (2.16) 

ijl 

aijl (p 2 r +j+i (m 2 c y in (4) * (r V - ™d . (2.17) 

ijl \ m bJ 

Terms with odd i, on the other hand, do not contain a logarithm. Note that we have 2m + 
2j + 21 — n + 1, which can be satisfied at any odd n. Thus we arrive at 

J d(v ■ p)^(v ■ p) 2 - p 2 P n {{v ■ p),p 2 )5 (n \p 2 - m 2 c ) (2.18) 
^Y\C im (p 2 Y n+3 - 2l ^ 2 (m 2 c ) l ]n (41 6(n) (P 2 ~ m ") ■ ( 2 - 19 ) 

ijl 

Hence IR sensitive terms can appear starting at n = 3 with the logarithmic dependence on m c 
of the Darwin term. For n > 3 we get from these terms a tree level contribution to the total 
rate of the form 

1/1 \ (rt ~ 3)/2 

r„ oc — ij — with n = 5,7,9... , (2.20) 

m b \ m cJ 

i.e., even a powerlike singularity for m c — >0. 

2.2 Explicit expressions for n = through 4 in the total width 

For n = the leading term reads 

P^ = \m 2 b {vp) (2.21) 
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which is an odd power of v ■ p and hence does not lead to IR sensitive terms as just explained. 
One should note that the subleading contribution to the leptonic tensor L^ v yields - upon 
contraction with the partonic hadronic tensor - a (subleading) term of the form 

P sub = -2m b (v ■ pf - p 2 m b ; (2.22) 

it leads to a contribution of the form m^ln(m^) in the phase space factor of the partonic total 
rate. We will focus eventually on those novel terms oc \jm k c that arise in leading order in l/m&; 
to low orders in A QCD it is easy to keep as well the terms stemming from the subleading pieces 
in L^. 

The next term with n — 1 is given by 

Pi = ^ G x ~^ (5p 4 + 7m 2 c p 2 - 20(v ■ pfp 2 - 10m 2 c (v ■ p) 2 ) , (2.23) 
where the hadronic tensor is contracted with the subsubleading part of the leptonic tensor 



" subsub 



. This yields again an ln(m^) terms upon integration over the phase space. 



For n = 2 the leading term of the leptonic tensor again contains only odd powers of (y ■ p) 
which do not generate any logarithms. 

At n = 3 the IR sensitive contribution is the Darwin term. Explicitly, we have 

{B(v)\b v {iD a ){iD y )(iD p )b v \B{v)) = ^M B p 3 D (g a p - v a v p )v^i> + 1) + • • • (2.24) 

from which we obtain 

P 3 ° ar = {3(p 2 -m 2 c f + 8(v ■ pf - 8p 2 (v ■ pf) . (2.25) 

Upon integration over (v ■ p) we arrive at terms with three types of prefactors: 

M 3 ln(g), <(P^(Q, < 2 26 > 

They have three derivatives with respect to p 2 from the (^-'-function, and hence the first term 
yields a ln(m^) factor which is the first infrared sensitive contribution. In the other two terms 
explicit factors of ml kill the infrared singularity and they thus remain finite for m c — >0. It is 
straightforward to check that in this way we end up with the correct prefactor for the infrared 
log in the Darwin contribution. 

The terms with n = 4 create again only odd powers of (v ■ p) and would yield infrared 
singularity for the l/m;,-subleading piece. Finally at n = 5 the following nine structures arise 

P 5 oc (v pf , (v ■ pfp 2 , (v ■ pfrrf , (v ■ pf{p 2 f , (v ■ pfm\ , (2.27) 

/ \2 2 2 / 2\3 / 2\2 2 2 4 

(vp) p m c , (p), {p ) m c , p m c 
Upon integration over (v ■ p) we obtain terms of the form 

« 4 K§)' (p2)3m HS)' ^ )2m HS)' p2m HS) (2 - 28) 

coming with five derivatives of the 5- function. All thus yield contributions of order 1/m 2 in 
the total rate. They will be addressed in detail in the next section. A similar consideration 
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extends in a straightforward way to higher orders where n > 5 emerge. These would generate 
terms inversely proportional to even larger powers of inverse charm mass. 

As a first resume we state that IR sensitive contributions - i.e., those singular for m c — > 
- unequivocally arise from the lower end of the integration over v ■ p (i.e. p — > 0) due to the 
presence of the non-analytic factor a/ (v ■ p) 2 —p 2 — \p\ in the integrand. For the dimension-six 
Darwin term they are of the form In (m 2 /m 2 ). Higher-dimension contributions exhibit an even 
stronger singularity, viz. powers of l/m c . Without extra gluon loops this happens first for 
dimension eight. 

The discussion above makes it clear that such IC effects emerge for the fully integrated 
width as well as for higher moments of the distributions. However the strength of the infrared 
singularity and the order in the heavy quark expansion it first emerges generally depends on 
the kinematic observable. 



3 nib ^> m c ^> Aqcd : charm as a dynamical quark 

In this section we shall consider an alternative way to describe IC effects. We now choose to 
integrate out the 'heavy' degrees of freedom only above the scale m c . This leaves the charm 
quark as a dynamical entity, much in the same way as would be required for light quarks in 
QCD, e.g. in b^uiu. The main difference that emerges is that, for inclusive decays, we have 
now to include four-quark operators explicitly containing charm quark fields. 

The charm quarks in b decay can act both as a hard and a soft degree of freedom. The hard 
component is treated the same way as described in the previous section. The matrix elements 
of the four-quark operators contain the "soft" part of the still dynamical quarks, which for now 
we treat nonperturbatively. Accordingly, we write the original QCD product of currents as 

(B\b{x)T u c{x) c(0)r„6(0)|B) (3.1) 
= (B\b(x)T u (c(x)c(O)) r p b{0)\B)>p + (B\b(x)T u c(x) c(0)T p b(0)\B) <fl . 

The product c(x)c(0) can be viewed as the Green function of the charm quark inside the external 
gluon field in the B meson, averaged over the field configurations present in the meson. This is 
an exact relation as long as the beauty meson has no charm flavor: it is a consequence of the 
Gaussian integration over the quark fields. This applies to the l.h.s. as well as to each of the 
two terms on the r.h.s. The decomposition on the r.h.s. merely reflects the different treatment 
used to describe these terms. 

The first term corresponds to the 'perturbative' (in l/m c ) calculation of the previous section. 
The second term has to be added now, since the charm quark is still a dynamic quark controlled 
by nonperturbative dynamics. The role of the intermediate scale \x is to draw the demarcation 
between the two dynamical regimes. 



Even though the first term in Eq. (3.1) is evaluated in the 'direct' way detailed in the 
previous section, the result differs due to the introduction of the cutoff \i. The precise form 
of how \i enters depends on the concrete way the (hard) separation is implemented. Let us 
mention that for tree-level calculations without extra perturbative loops it is sufficient and 
convenient to simply integrate the distribution with the constraint 

< q 2 < {m b -fi) 2 . (3.2) 



The upper bound above effectively introduces the separation scale in Eq. (3.1) when the corre- 



lator is integrated over the phase space to obtain the inclusive width or its moments. 
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Evaluating the first term in Eq. (3.1 ) in the previous section gave contributions IR sensitive 
to the charm mass. Taking the formal limit m c — > separately term by term in the expansion 
would have yielded divergent expressions. With a nonzero p ^> Aqcd this changes: all individual 
terms remain regular by themselves at m c — > 0. For the role of the infrared regulator is now 
taken over by the Wilsonian cutoff. 

This is evident on general grounds: the terms IR-singular for m c — > came only from the 
soft charm configuration with momentum p < m c - the domain now excluded. Alternatively, 
this can be traced explicitly in the formalism of Sect. 2. For instance, with the constraint of 
Eq. (3.2) the lower limit of integration in (v ■ p) rises from ^JpP to p +2 ^^~ M — A*; therefore 
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III 



log -fit m Eq. (2.8) turns into \og-Jt. All the integrals become analytic functions of m c at 



m c <^p] the logs and inverse powers of m c are replaced by those of the cutoff mass p. 



The second term in the r.h.s. of Eq. (3.1) has a smooth m c — > limit. Although it may 
have soft "chiral" singularities when both m c and one of the light quarks become massless, 
the expectation values should remain finite; only higher derivatives with respect to the charm 
mass may have singularities if m u or vanish. We note that this expectation value, being 
regularized in the ultraviolet, is well-defined and these conclusions hold with no reservations. 

We may briefly discuss at this point the parametric dependence of the total semileptonic 
width r s i(6 — > c) on m c when the latter is small. As long as m c ^> Aqcd holds, we have 

A 3 2 A 3 / A 2 \ ^ 

non-analytic nonperturbative terms at order 1/mf scaling like In ^| , f ^2° J with 

k>0 (in general, odd powers of m c also emerge). Each of these terms separately are singular 
at m c — > 0. The leading term is driven by the Darwin expectation value; it represents an 
IR singularity in m&/m c which, in principle, is observable, at least at large mj and sufficient 
accuracy. 

The expansion in Aqcd / Tn c makes, however, sense only as long as charm remains heavy on 
the scale of QCD dynamics. At lower m c the successive terms with higher k would formally 
dominate. The whole function of m c stabilizes at m c < Aqcd and approaches a finite value at 
m c ^0. A model for such a behavior can be given, for instance, by 

(here A is a strong interaction mass scale parameter of the order of Aqcd) which, expanded in 
1/m^, would yield the whole series in l/m 2 c . The actual coefficients for the l/m 2k terms may, 
of course, be different, and they can be calculated in the OPE along either road. 

Returning to the OPE analysis proper, we can bridge the two ways of accounting for the 
nonperturbative charm effects by looking at the //-dependence of both terms in Eq. ( 3.1[ ) (more 



accurately, when it is integrated to obtain the inclusive probability). Their sum must be /x 
independent which provides useful relations. The following note should be kept in mind. 

The form of the /i-dependence is determined by the regularization scheme. In the Wilsonian 
procedure with a hard cutoff the leading log dependence is accompanied by power terms. The 
resulting dependence is qualitatively different for /i<^m c and for /i>m c . 
(i) When p, is taken small compared to m c , it enters only as a small power correction, ~ (p/m c ) 1 ; 
in the formal limit p—>0 (m c ^> Aqcd fixed) the last term would vanish, the whole correction to 
the width is then given by the first term. Raising p up to m c and above moves the correction to 
the width from the first term to the second. The calculation of Sect. 2 represents the evaluation 
of the last term contribution in the Aqcd /m c -expansion for p~^>m c . Such a calculation makes 
sense only as long as charm is sufficiently heavy. 
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(ii) At fi^>m c the situation is different. Beyond the leading term On In ^ the /x dependence 
is suppressed by powers of ml/fi 2 . In the second term in Eq. (3.1) these appear as a residual 



dependence of the high-dimension terms on the ultraviolet cutoff for integrals intrinsically 
convergent at the scale m c . In the first term it shows up as now small coefficient functions of 
higher-dimension operators - which would normally be saturated at soft charm configurations 
p ~ m c - proportional to l/m^l//i 2fc , instead of 1/mf without a cutoff. It is then 

convenient to assume /i>m c and neglect these terms altogether. 

A short comment is in order on how the separation would look like in the often adopted 
dimensional regularization scheme (DR) where no powerlike dependence on renormalization 
scale ever arises. For small m c the only /i-dependence within DR enters through ln/i in the 
Darwin operator. In such a scheme, however, charm quarks must be treated as massless; the 
dynamic four-quark operators have to be renormalized in the UV likewise in the way of DR. 
The requirement to set m c = then is rather clear. For keeping m c finite would destroy naive 
DR in a direct calculation: it yields a finite result for the first term around D = 4 because the 
potential IR singularity at D = 4 is regularized by a non-zero charm mass. As a result, with 



m c y^0 calculating the first term in Eq. (3.1 ) with DR precisely reproduces the total contribution 



of Sect. 2. Consequently, in DR at m c ^ the last term in Eq. (3.1) should be considered as 
vanishing regardless of the concrete value of the charm mass; this is not a very physical result, 
at least for a relatively light quark. 

This is a rather typical feature of dimensional regularization. DR may provide a convenient 
highly efficient technical tool to analyze the case of massless final state quark or when a non- 
zero charm mass can be neglected. However, it would require additional matching procedure, 
usually order by order in m c if the charm mass dependence is important. 

In logarithmic terms - as for the Darwin operator in the integrated width - the charm 
mass plays the role of a renormalization point. Therefore, as discussed in Ref. [2], the infrared 
dependence on m c in the 'conventional' calculation along the 'first' road must match the UV 
dependence of the corresponding four-quark expectation value given by p 3 D In A similarly 
dual description applies also for the terms scaling like inverse powers of m c . They can be 
calculated conventionally assuming charm to be heavy following the route of Sect. 2. Alter- 
natively, they can be obtained as the corresponding pieces of the four-quark expectation value 
(normalized at /i>m c , Aqcd)- The results over either road must be identical whenever one is 
in the domain where the expansion can be applied. The first road (without implementing a 
cutoff) is, of course, justified only for m c 3> Aq CD . The second route is formally valid for an 
arbitrary hierarchy between m c and Aqcd _ however, we do not have the means to calculate 
the expectation value through gluon operators without charm fields when charm becomes light. 

In order to make the above mentioned correspondence explicit, we will address the 1 / mf l/m™ 
terms where at tree level only even n emerge. To this end, we consider the contribution of the 



second term in Eq. (3.1) involving the explicit charm quark operators. Inserting it into the 



hadronic tensor we get 

2M B W^ = J d A x exp(ip • x)(B(v)\b v (x)T u c(x) c{Q)^ M^)\ B ^))\^ (3-4) 
= J d 4 x exp(i P 'x)(B(v)\b v (0)r v c(0) c(0)TM^)\ B ( v ))U 
+ J d 4 x ex P (ip ■ x)x a (B(v)\[d a (b v (0)T ,0(0))] c(0)r it b v {°)\ B ( v ))\» + ' ' ' 
where the ellipses denote operators of dimension eight and higher. Extra derivative for (renor- 
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malized) expectation values can bring in a factor of m c or \i at most; it can be traced that 
powers of x a translate into powers of l/m&. Therefore, the higher terms in this expansion 
yield higher powers in the l/m& expansion. We then can focus on the first term which is the 
conventional D = 6 four-quark operator. The corresponding hadronic tensor reads 



2M B W^) = (2n)H\p)(B\b(0)rA0)cmA0)\B)» 

+ (2vr) 4 (-i^-) S\p) (B\d a bT u c(0) c(0)r M &(0)|a)„ + ... (3.5) 



9p a, 

and we retain only the first term. Note that the 5-function projects out the leading term of the 



leptonic tensor (2.10). Furthermore, this contribution is localized at p = - hence at p 2 = 
and v-p = 0, in agreement with the findings of the last section. As a consequence, step- functions 



in Eq. (2.1) become superfluous. The last relation completing the arithmetic part is 

/d(v-p) dp 2 a/ (v-p) 2 —p 2 S 4 (p) = — . (3.6) 
27T 



del 

Omitting QCD corrections, the relevant diagrams for calculating W^\ u in Eq. (3.5) are the 
one loop diagrams involving the charm-quark loop with an arbitrary number of external gluons. 
These diagrams have been considered already in [lj . It is advantageous to first perform a Fierz 
rearrangement of the four quark operator according to 

bY v c cT% = -h a T p b p cpT a c a [-ie^° + <f V" + <T<T - Q^sT] , (3.7) 

where a, (3 are color indices, and the minus sign comes from anticommutativity of the quark 
fermion fields. 

We may construct the charm mass expansion of this charm loop in the external gauge 
field and take the average over the S-meson state. There are a few subtleties related to this 
procedure, since the two charm quark operators are taken at coinciding space-time points. As 
in pQ we may start from the time-ordered product of two charm quark operators at displaced 
points, which amounts to consider the conventional charm propagator in an external field. 
This propagator is generally gauge dependent, yet in the end this is compensated by the same 
displacement in the 6-quark fields. For constructing the short-distance expansion of the Green 
function the fixed-point gauge is convenient. 

The limit of coinciding points in the charm Green function is formally divergent, yet gauge 
independent. Subtracting the free Green function (this piece is accounted for in the purely 
partonic width) we end up with a mild log divergence present in the vector current, proportional 

to [D^G^y. 

2 1 /A? 



(c a Yc,) A = I j^y 2 In ^ j [D K , G KV \p a + ■ ■■ (3.8) 

where G pu is the gauge field strength tensor; A uv should be identified with /i in this context, 
and the ellipses denote finite terms to be considered below. As discussed in |2j , this ultraviolet- 
singular log matches onto the infrared piece of the conventionally calculated Darwin coefficient 
function. The contributions from the axial-vector cc current are convergent ab initio. Examples 
of the Feynman diagrams are shown in Fig. [Tj 
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b b b b 



Figure 1: Diagrams illustrating calculation of the charm loop in the external field. The wavy 
lines generically reflect insertions of the external gluon field. 



The lowest finite terms yield 1/mf 1/m 2 . contributions given by 



4&n 2 m 2 c 



2\[D k ,G kX ],G vX } + 



D K ,G 



G 



Pa 



+ ■ 



2407r 2 m 2 



(13 [D K , [G Xu ,G XK ]]+8z [D k , [D x ,[D x ,G ku ]]] 
-Ai [D\[D K ,[D x ,G h 



(3.9) 
(3.10) 



This assumes /i>m c and neglecting power terms ~ [m c j \x) . Inserting this into (3.7) we end 



up with dimension-eight b... b operators with gluon fields; their coefficient functions compared 
to the partonic D = 3 operator bb are proportional to l/mf 1/m 2 ,. 

A closer look into the calculations of the charm loop in the external field reveals that the 
resulting expressions exactly parallel those in Sect. 2 once the leading-order approximation 



Eq. (2.10) is adopted and only non-analytic terms according to Eq. (2.7) are retained. This 



occurs before the full integration is performed, when one takes the integral over the timelike 



component of the loop momentum by the residues at p 2 = m? c) for each power term in the 
expanded propagator. 

In fact, the technique of the loop calculation in the external field itself allows to derive 
a number of relations which strongly constrain the form of the operators which can appear 
in such an expansion. It has been presented in detail in the first part of Ref. [6]. These 
relations ensure that the result has always the form of multiple commutators. This sharpens 
one observation made already in Sec. |2j There we noted that the partonic matrix elements of 
IC contributions necessarily have to be at least one-gluon matrix elements. From the arguments 
given in this section we conclude that intrinsic charm contributions involve only gluon fields 
and their derivatives; there will be no derivatives acting on the beauty quark fields that would 
generate a dependence on the 'residual momentum' of the decaying 6-quark. 

We have verified though explicit calculations of the l/m\ 1/m 2 terms that these two ways 
to calculate l/m c -singular terms yield the same result for the operator expansion. 

At order ljm\ljm 2 c the result can thus be expressed through five operators, which are 
determined by two contributions to the axial vector current and three contributions to the 
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vector current: 



2M B U 



2M B ~h = (B\b v \iD K} [iD x , [iD\iG^]]\b v \B) v a 

2M B f 2 = (B\b v [iD x , [iD K , [iD\iG Ka ]]\b v \B) v a 

2M B f 3 = (B\b v [iD K , [zG Xa ,tG XK ]]b v \B) v a 

(B\b v { [iD p , iG pX ] , iG Sl }{-ia a p)b v \B) 



x 



(g Xa g 



2M B ~h = (B\b v {[iDe,iG aX ],iG n }(-icr a p)b v \B) 



(3.11) 
(3.12) 
(3.13) 

(3.14) 

(3.15) 



The contributions originating from the axial current yield spin-triplet operators, while those of 
the vector current yield spin singlet operators. 



4 Quantitative estimates of Intrinsic Charm 

To estimate the IC contributions we follow the lines of [T] and apply the "ground-state factor- 
ization" approximation to the matrix elements. 

The effect on the total rate has been considered already in [I], yielding a reasonably small 



contribution. In fact, the total rate, expressed in terms of the operators given in (3.11) reads 
as 

3 2 

i ~~ 



m 3 b m 2 c 



To 



-8/i+4/ 2 -13/ £ 



1 2 



(-2/4 - h 



(4.1) 



2 15 V J± -V 2 3 

Following the way to evaluate the expectation values suggested in Ref. [Tj we obtain numericalljj^] 
,A^0.31GeV 5 , / 2 ^0.25GeV 5 , / 3 -0.14 GeV 5 , / 4 ~0.34 GeV 5 , / 5 ~-0.40GeV 5 (4.2) 



which lead toE] 



ST 



(0.7%) x r 



Parton 



(4.3) 



This numerical estimate should not be considered bullet-proof. While the individual matrix 
elements are predicted with reasonable confidence in their signs and magnitudes, we are faced 
with a set of terms with different signs. Thus cancellations will in general occur among them. 
Their degree may depend on the numerical accuracy of the applied ground-state factorization, as 
well as on the precise values of the lower-dimension expectation values fi 2 , p z D and p LS . A more 
elaborate discussion of the nonfactorizable effects will be presented in a separate publication. 

Finally we note that for the moments the situation is different. Contributions that introduce 
an IR sensitivity to the charm quark mass in the total rate may become regular for the moments. 
This becomes evident if we consider moments of the partonic invariant mass such as ((p 2 —rn 2 c ) n ) 
which remain regular as m c — > till higher orders in Aqcd- 

2 We have found a discrepancy with Ref. p] in the overall factor for one of the expectation values. It does 
not produce a noticeable numerical change for the correction to the width, however. 

3 Ref. [lj included an additional phase-space suppression factor for the IC kinematic of (l—m c /mb) 2 - Based 
on the operator analysis we can show that actually it is absent in the case at hand. 
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Having at hand the numerical estimates of the higher-dimension expectation values allows us 
to refine the model (3.3) for the charm-mass dependence in the IR regime. Since the l/(m\m 2 c ) 



corrections calculated in this ansatz are fixed in terms of parameter A, we assign the latter the 
value which would reproduce these leading corrections. This yields 



A 2 = M 2 
and the model predicts 



h + 13/s 



h 



A 

5p% 10p 3 D ' 40p 3 D 12p 3 D 24pl 



(0.7 GeV) 2 , 



(4.4) 



( _^ +w V ) _ {B |5 7a(1 _ 7b)cgT/j(1 _ 7b)6 | B)/ 



P 3 n 



In 



4tt 2 "m^ + Mf 

see Fig. [2J In this case the correction to the width at m^Cmj takes the form 



(4.5) 



5r s 



mi 



In 



mi 



77 



m 



l + Ml 48 



(4.6) 



where the constant term accounts for the explicit UV contribution in this limit. This model 
illustrates to which extent the charm quark may be considered heavy in this context. 



(bccb) I GeV 3 



-0.005 - 



-0.010 - 



-0.015 



-0.020 - 




rn c I GeV 



Figure 2: The WA expectation value 2M~{B\b^f (1— 75)007 (1— 7 5 )6|_B) M as a function of charm 
mass: the leading contribution logarithmic in m c (blue dashed), including additionally the l/m 2 



power correction (orange dotted), and the complete behavior according to ansatz Eq. (4.5) (red 
solid). We assume the value p\ = 0.15 GeV 3 and the normalization point [i is taken 4.6 GeV. 



It should be acknowledged that Eqs. (3.3), (4.4) and (4.6) are only a reasonable model for 
the effects at intermediate to small m c . In fact, the size of the effective mass scale M* as 
determined by matching the leading charm power corrections may not be fully universal: it 
depends on the particular Lorentz structure of the weak vertices. It would likewise differ, say, 
in weak B* decays - although numerically such a variation in M* would be insignificant. 



5 Facing up to Weak Annihilation 

Our preceding discussion can suggest a dual description of IR charm effects that allows insights 
into a higher-order OPE calculation of the heavy-to-light case B — > X u £i?£. More specifically the 
limiting case m c — >0 is of relevance when treating the heavy-to-light case beyond order l/m 2 . 
The inherent IR divergence of the standard calculation for m c /mfo— >0 underlies the importance 
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of the four quark operator matrix element in the heavy-to-light case, which is usually called WA 
contribution. More precisely, here it corresponds to its valence-quark insensitive piece which 
affects semileptonic B + and B° decays equally. We refer to it as non-valence WA. 

The important question here is the potential scale of the non-valence WA contributions, 
once the effect of the Darwin operator has been separated outj^] A direct computation based 
on the method discussed above is not possible. Yet one may try to estimate the natural scale 
expected for WA by approaching the massless case from the heavy quark side. For this purpose 
we use the model suggested in the previous sections. In spite of its admitted oversimplification, 
we can be sure that the difference between this ansatz and the actual QCD contribution remains 
finite at any m c including the limit m c —>-0. The model has the advantage of being sufficiently 
accurate when extrapolating down from the side of intermediately heavy charm quarks. 

It is therefore plausible that setting m c = and applying this model we do not stray far away 
from the leading effect of the non- valence component of WA in b -^ulv decays (which is its 
only contribution in B d decays) Q this guesstimate, though, cannot be validated in the context 
of the l/m c expansion examined here. In some sense such an assumption implies that no 'phase 
transition' occurs when going down in mass from heavy to light quarks. We know such a phase 
transition takes place in the QCD vacuum, yet it may not necessarily be important for the 
expectation values over B meson states. One may a priori expect larger WA effects in b^ulv 
coming from its flavor-specific piece manifesting itself in the decays of charged 5-mesons. 

Using the model of the previous sections to interpolate between the regimes of heavy and 
light charm we get an estimate 

1 -{B d \bj(l- l5 )uuj(l- l5 )b\B d ) ~ -£*L\n( « -0.005GeV 3 , (5.1) 



with fi denoting the normalization scale. This is an educated guess and cannot guarantee to 
yield even the correct sign of the effect at // < 1 GeV. The negative sign physically means that 
the propagation of the soft u quark (projected onto the spin state specified by the Lorentz 
structure in question) is suppressed compared to free propagation. Taken at face value, the 



WA expectation value in Eq. (5.1) would yield the isoscalar shift in the semileptonic B decay 
width 

6TT(^«)_. ^4^ + 1 U_ . 015 . (5 . 2) 



T sl (b^u) ml 

A more refined way to assess non- valence WA for b — > u iv would be to consider the real 
massless case for the final-state quark, m u = 0, yet to introduce an IR cutoff via the kinematic 



restriction Eq. (3.2). The first term in Eq. (3.1) representing the £ UV piece of the buub 



expectation value from the domain of quark momenta above \i is then calculated in the direct 
way of Sect. 2 (with m c — > 0). The result is expressed in terms of the same five expectation 
values, with the coefficients scaling as l//i 2 ; yet they would combine to yield in general a 



4 The difference due to the valence part of WA in b^uiv can be experimentally probed by analyzing the 
difference in the semileptonic spectra spectra of charged and neutral B decays. 

5 Paper U which first analyzed the effects of generalized WA in semileptonic decays, focussed on the differ- 
ences between mesons with different spectators and therefore explicitly subtracted the bu ub expectation value 
in Bd from that in B~ . Following an earlier classification of the preasymptotic power corrections to the inclusive 
widths it referred to this as the spectator-dependent correction, a terminology continued to a number of later 
publications. The valence and non-valence effects separately were considered in Ref. [7] where valence effects 
were sometimes also called 'spectator' contributions. The importance of the non-valence WA for light quarks, 
although conjectured already back in 1994, was emphasized in Ref. [8 . There the 'singlet' and the WA proper 
pieces referred to the average and the difference of the valence and non- valence components. 
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different combination of the operators and, consequently, a different number. Evaluating the 
result with /i^0.6GeV would provide an estimate of the minimal natural scale of non- valence 
WA. 

Even such an estimate would be admittedly incomplete; beyond its lack of precision in 
evolving the \i dependence to as low a value as 0.6 GeV, the total WA should also include 



the (-B|6MM6|-B)<o.6GeV ; the last term in Eq. (3.1). The contribution from the low momenta 



plausibly exceeds the numerical estimate above, and may even change the overall sign. However, 
it would be unnatural to allow the contributions from physically distinct domains of low and 
high momenta to show significant cancellations. Therefore, we would view the thus obtained 
estimate as a firmer lower bound on the scale of non- valence WA in b^uiv. 

One may anticipate a potentially more significant effect for the 'valence' part of WA which 
describes this effect in T(B + ~^>X u tv). On physical grounds we expect this contribution to 
contain a piece independent of the non-valence WA and not related to something that can be 
traced from the b^-ciu decays in the limit of small charm mass. In the formal derivation of 



relating c(x)c(0) in Eq. (3.1 ) to charm Green functions it would be associated with an additional 
term. That term appears through the contraction of the quark fields with those of the same 
flavor whose presence is required in the interpolating currents to produce the initial and to 
annihilate the final B meson state. 

On the other hand, as pointed out in Refs. [U [9] and exploited in later papers [8], [7], one 
may infer certain information about WA, and in particular about valence WA from the D me- 
son decays. The charm quark is marginally heavy to apply precision heavy quark expansions 
to its decays. Extrapolation from charm to beauty may thus be semi-quantitative at best, yet 
it should provide some constraint on the expected scale of WA's physical implementation. In 
particular, we want to point out that the preliminary CLEO-c [10] data on the semileptonic 
branching fraction of D s indicate the presence of a destructive spectator-related WA contribu- 
tion of around 20%; it must be the result of non-factorizable four-quark expectation values. In 
interpreting this observation one has to address the question of 5'L r (3)-breaking in the lead- 
ing nonperturbative corrections described by the kinetic and chromomagnetic operators (and, 
possibly, in the Darwin expectation value). 

While semileptonic 'valence' WA cannot occur at all in D° decays, it can contribute in D + 
and D s decays on the Cabibbo suppressed and Cabibbo allowed levels, respectively. Assuming 
possible SU(3) breaking to be under control in WA proper, we conclude that the observed 
difference of the total semileptonic widths for D s and D° 

r sl (D°) ~ Br sl (D°) r Ds ~ ^ ^ 

must be dominated by valence WA in D s . To describe the pattern on the WA effects in D 
mesons we may adopt the nomenclature for the generalized 'annihilation' correction following 
Ref. jjj: The valence WA for a particular transition c^qiu (q = s or d), WA™ 1 refers to the 
difference in the matrix elements between D q and D°. The nonvalence WA^ val is directly the 
expectation value in the (non-strange) state D° containing no valence d quark. To allow for 
SU(3) asymmetry we also have to distinguish WA^ val from WA^ val(s ' ) where it is considered 
in the strange D s state. The WA operators above may be general products like cqqc, either 
local or nonlocal. (For the decaying quark being heavy enough, like in 5-mesons, it would be 
sufficient to include only the leading local four-quark operators.) 
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With this convention, we in general have in D mesons 
r s i(-D + ) — r s i(D°) = sin 2 fl c -WAf 



T sl (D s ) - T sl (D») = cos 2 9 C • WAJ al - sin 2 9 C [WA^ al - WA° val(s) j + A su{3) . (5.4) 

By introducing the subscript marking the d or s flavor we have explicitly allowed for the SU(3) 
breaking in the expectation values due to the different spectator in a meson or in the light quark 
field flavor in the corresponding operator. We have still neglected the explicit short-distance 
SU (3)-breaking oc m 2 s in the coefficient functions emerging due to the larger m s in the hard 
quark Green functions; it is expected to be strongly suppressed. A$u(3) in Eq. (5.4) therefore 



denotes only the shift related to the 577(3) violation in the (flavor-singlet) nonperturbative 
expectation values between the strange and non-strange heavy meson states. The analysis 
suggests that these effects are numerically suppressed for the kinetic and the chromomagnetic 
operators and should not exceed 5% level in the widths. Then the bulk of the difference in 



Eq. (5.3) should be equated with the valence component of WA, at least if SU(3) violation in 
it is not too strong. 

Translating these relations for WA from charm to beauty is associated with significant 
uncertainties due to a potentially poor representation of the contributions to the inclusive 
width for charm by the (truncated) OPE; for m c is manifestly not large enough for a precision 
treatment. Relating WA for B decays to the expectation values of the D = 6 operators can be 



done with acceptable theoretical accuracy. Yet expressing the WA contributions in Eqs. (5.4) 
through the analogous local expectation values in D mesons is subjected to large corrections 
- first of all from the corresponding higher-dimension operators with additional derivatives. 
Related to this is the short-distance 'hybrid' [TTl[T2] renormalization of the operators in question 
from the scale of charm to beauty; in this case it may be even not fully perturbative. 

Bearing in mind these potential caveats we nevertheless use this line of reasoning to estimate 
the expected size of the valence WA in the semileptonic b— >u width of charged B meson: 



r(B^x ]ight eu) 

which is similar in magnitude, yet still below our estimates for the minimal scale of the non- 
valence WA. 



6 Summary and Outlook 

The main result of this study is that the OPE for inclusive B — ► X c ivi contains terms with 
an infrared sensitivity to the charm-quark mass. Although this has been known, a complete 
discussion of these so-called "intrinsic charm" contributions had not been presented. We have 
given such a discussion here in the context of two theoretical frameworks or 'roads' for removing 
c quarks from the dynamical degrees of freedom that a priori appear different, yet in the end 
yield identical results. 

We have shown that starting at l/mf the standard OPE for B— >X c li>i exhibits terms of the 
form l/m™x 1/m™ where at tree level only even n and odd m appear. The matrix elements of 
local operators parametrizing their nonperturbative input always contain gluon-field-strength 
operators and their covariant derivatives; in turn the residual momentum of the b quark does 
not enter here. 
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We have performed a detailed analysis of the contributions of the form 1/mf xl/m? c at tree 
level, which is needed to complete the OPE calculation of B — > X c ti>i up to order l/mf, since 
parametrically 1 jm\ x Aqcd/^c i s °f the same order. The numerical estimates confirm the 
results presented in Ref. []]. 

The conclusion for B —> X c lvi is that a calculation to order has to include also the 

terms of order x l/(m 2 c ) k \ furthermore, including radiative corrections one obtains also 

contributions of the order 1/m™ x a s (m c )/m k where k can also be odd. The lowest terms of 
this kind are of order l/mf x a s (m c )/m c and have been considered in Ref. pp. 

These effects are of considerable theoretical interest with respect to subtleties that can arise 
in nonperturbative dynamics, yet they go beyond it towards more pragmatic goals: they help 
to validate the goal of reducing the theoretical uncertainty in extracting \V c b\ from B—>X c £i> 
to the 1% level; achieving such a goal is of obvious interest for the theoretical treatment of B 
decays - yet also for a proper interpretation of the ultra-rare decays K — > tt vv. Their amplitudes 
have been calculated with high accuracy in terms of m c and V t *V t d |13j . Their widths thus scale 
with Vjh, and the error on the latter is at least a large component in the stated overall 2% 
uncertainty. 

In b^ulv decays the straightforward calculation of the higher-order power corrections to 
the total width beyond order 1 jm\ would yield terms which diverge power-like in the infrared. 
Yet this does not mean that one cannot go beyond \jm\ order here. The analysis shows that 
to calculate r sl (6— >) without extra a s -corrections it is sufficient to introduce the corresponding 
WA four-quark operators, and then one should simply discard all the terms formally having 
inverse powers of m u . 

We have found in the process that analysis of IC effects can inform and focus our thinking 
about the possible impact of WA in the heavy-to-light transitions B — > X u Iv - and the ex- 
traction of \V u f)\ there - and on its relation to charm decays. The IC effects for the inclusive 
distributions are conceptually similar to (generalized) WA corrections extensively discussed in 
connection to the lifetimes of heavy flavors since the late 1970s. Since the usual beauty hadrons 
do not contain valence charm quarks, we deal here with the case of non-valence WA contribu- 
tions first noted in Ref. [3]. A profound difference with the conventional WA for light quarks is 
that charm quarks, even soft ones, to the leading approximation can be viewed perturbatively, 
and the nontrivial strong dynamics affect its propagation at the level of power corrections 
(Aqcd / 'm c ) k , while for light quarks this would not represent a parametric suppression. Never- 
theless, approaching the case of conventional WA with light flavors from the heavy-mass side 
and using a model which naturally interpolates between the regimes of heavy and light quark 
we get an estimate 

* ./ - 1 « -0.015 . 6.1 

This result should be viewed as an educated guess rather than a real evaluation; one cannot 
count even on the firm prediction of the sign. It also leaves out the more intuitive valence 
WA which, as a matter of fact, historically gave the phenomenon its name. (Yet, as clarified 
in Ref. [12J, its interpretation in the presence of strong interactions is more subtle and may 
include interference-type contributions which allow for the net correction to the width even to 
become negative.) 

We have used the recently reported [TU] measurements of the D s semileptonic fraction to 
estimate the significance of spectator-related WA in the KM-suppressed semileptonic width of 
B + . Taken literally, the correction turns out to be close to the non-valence case Eq. (6.1) or 



16 



even somewhat smaller, but still destructive: 



-(0.005 to 0.01) . (6.2) 



Since these contributions populate the kinematic domain of small hadronic invariant mass 
and energy of the final hadron state, they could have an amplified impact on the existing 
determinations of V u b from B — > X u £p£. 
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